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INTRODUCTION 


During  the  past  two  decades  significant  <•  d vances  have  been  made  in 
the  theory  of  non-linear  elasticity.  One  of  the  potentially  most  useful 
areas  of  research  in  non-linear  elasticity,  from  an  engineering  viewpoint, 
has  been  the  development  of  small  deformation  equations  superposed  on  a 
state  of  finite  deformation  [1].  A  logical  extension  of  these  results 
would  be  to  adapt  these  equations,  by  means  of  an  averaging  procedure, 
to  describe  beam  and  plate  behavior.  For  example  the.  spirit  of  this 
suggestion  has  been  carried  out  for  plates  by  Herrmann  and  Armenakas  [2], 
who  employed  a  variational  method  to  derive  their  equations.  An  alternate 
procedure,  which  is  used  in  the  present  paper,  is  to  start  with  the  al¬ 
ready  developed  equations  of  non-linear  elasticity  and  then  to  sequentially 
perturb  the  equations,  integrate  thru  the  thickness,  and  finally  introduce 
displacement  simplifications  to  obtain  a  linearized  beam  theory  subjected 
to  an  initial  non-uniform  state  of  stress. f  The  advantages  of  this  proce¬ 
dure  are  two-fold.  Firstly  it  permits  some  physical  understanding  to  be 
attached  to  the  development  of  the  governing  equations  and  secondly  it 
provides  a  check  on  the  equations  (one-dimensionalized)  developed  by  the 
orevious  variational  technique  [2]  and  its  associated  assumptions.  As 
will  be  seen  the  two  sets  of  governing  equations  are  at  wide  variance, 
and  since  the  present  technique  is  easily  followed  and  checked  one  tends 
to  believe  the  present  results.  The  details  of  this  procedure  are  de¬ 
veloped  in  the  following  sections. 


t  Obviously,  this  same  technique  may  be  used  to  generate  a  plate  theory, 
which  will  be  presented  at  a  later  date. 


THE  NON-LJiEAR  EQUATIONS  OF  ELASTICITY 


Referring  to  t-ijure  1  and  following  accepted  procedures  [3],  the 
stress  vector  equation  of  equilibrium  of  the  deformed  body  is  given  by 


3  oi*/3  x.  +  X*  =  0 


v  a.*  '.i.  t*  e  stress  vector  referred  to  the  undeformed  i  face  area 

an*  '  is  rhe  hd.  ‘ir-ce  v rctor  referred  to  the  undeformed  volume.  Re- 
solv.'-g  a. .*  i"  aon-orthogona i  lattice  vector  (G.)- directions  yields 

•  J 

sf*  -  v  h 

where  a-.*  are  the  Trefftz  components  of  stress,  referred  to  the  unde- 
* 

formed  i^1  face  area,  which  can  be  shown  to  be  symmetric  [3].  Since 

OD  4- 

G.  =  r-  where  R  is  the  final  5  Late  position  vector  it  is  seen  that 
j  dXj 

-*■  3u  -*• 

G .  =  (6-  +  ■—)  1  (3) 

j  v  JS  3x.'  s  v  1 

J 

where  i  is  the  unit  vector  in  the  s^  orthogonal  cartesian  direction. 

•+■  *» 

Decomposing  X*  into  its  i  components 

X*  =  X*  i$  (4) 

and  putting  the  results  of  (2),  (3),  and  (4)  into  (1)  yields  the  follow¬ 
ing  scalar  equations  of  equilibrium. 


»  U  LI 

3x7  C  (6ds +  1  +  xr  *  0 


t  Note  that  R  =  r  +  u=(X$  +  u$)  i$  ,  therefore  =  (6.$  +  j^-)  i$ 

J  0 
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Referring  to  Figure  2  where  p*  is  the  prescribed  traction  referred  to  the 
undeformed  oblique  face,  and  noting  that  dA  (i.  •  n)  =  dA.j/2,  the  equili¬ 
brium  of  the  deformed  tetrahedron  is  given  by 


p*  =  (i  •  •  n)  a-*  =  n.  g- 


(6) 


and  by  using  (2)  and  (3)  and  defining  the  i  components  of  p*  as  p*  , 


(6)  becomes 


3u. 


»£  ■  V  "i  <5is  *  w:> 

J 


(7) 


Equations  (5)  and  (7)  are  the  desired  non-linear  equations  that  describe 
the  equilibrium  condition  and  the  surface  tractions  respectively.  If  the 
extensions  and  shears  are  small,  the  final  areas  and  volume  are  equal 
respectively  to  the  initial  areas  and  volume  so  that  a..*  =  a- .  ,  X*  =  X 


U 


and  p*  =  p  where  a. X  ,  and  p  are  the  actual  stresses,  body  forces, 

•>  l  J  ^  •> 


and  surface  tractions  respectively.  These  approxi motions  are  adopted  for 
the  work  that  follows,  so  that  the  equilibrium  equations  and  the  surface 
traction  equations  become, 

•v  8u_ 


sr  r  <6js  +  ^ 


■)  a. .  ]  +  X  =0 
ij  J  s 

(8) 

8u_ 

+  air> 

(9) 

THE  PERTURBED  EQUATIONS 


Following  a  technique  described  by  Bolotin  [4]  the  following 
quantities  are  introduced 


u  =  u  +  u 
s  s  s 


(10) 


(11) 


c. .  =  a • •  +  a-  • 
1J  1J  U 


Ps  =  Ps  +  A  ps  +  Ps 


X,  =  Xc  +  A  X$  +  Xs  -  p 


where,  for  example,  u$,  u$,  and  ug  represent  the  final  displacement,  the 


initial  displacement,  and  the  perturbation  displacement,  respectively. 


The  terms  A  p  and  A  X  represent  the  change  in  the  initial  tractions 


and  body  forces,  respectively,  due  to  the  perturbation  displacement. 


Since  the  final  equilibrium  equations  and  traction  expressions  have  forms 


similar  to  (8)  and  (9)  one  has 


_  C7  U 

air  [  <sjs  +  sir*  5ij  ]  +  xs  =  0 


n  =  a. .  n.  (6 .  +  7~) 
Hs  u  i  js  9x_. 


and  upon  entering  (10)  thru  (13)  into  (14)  and  (15),  and  using  the  results 


(8)  and  (9)  one  has 


"  (t,ij  axT*  +  air  £  “ij  (5sj  +  at +  at*  *  +  xs  +  4  xs  *  p  “5F  =  0  (16) 

1  W  1  J  J 


9u_  9u 


9u  _  9u  9u 

Ps  +  4  Ps  *  C  °ij  STT  +  ®1j  (5js  +  3x7  +  3x7*  *  ni  (17> 

\J  vl  J 


Linearizing  the  perturbation  quantities,  one  neglects  the  term  a^.  , 

0 

and  assuming  the  initial  displacement  gradients  to  be  small  [4]  one  neglects 
9u„ 


the  term  r-=-  .  This  results  in  the  following  set  of  relations. 

oX . 


a  9uc  a  _  -  32uc 

sir  (°fo  a?:1  +  at  <0is>  +  XS  +  4  XS  -  p  IF  =  0 

I  J  * 


Ps  +  A  ps  -  [  aij  9x -  +  ais  ^  ni 

J 


THE  BEAM  PERTURBATION  EQUATIONS 

Assuming  a  two-dimensional  problem  (plane  stress)  in  the  xi  -  X3 
plane  (18)  and  (19)  become  upon  expansion* 

-i_  (a  IHl)  +  _L_  (fj  Z&L)  +  _L_  la  \  +  3  (a  SlTi  N  8an  9a31 

3Xj  1011  3xi 1  3xi  101  3  3x3  3x3  {031  3Xl}  3^T  (033  377}  +  Z7T  +  3^~ 

+  7,  +  A  X,  -  p  ^  =  0  (20) 


A  <°»  !£•>  +  sir  (°>3  !&>  +  5TT  (®»1  4£>  +  A  («»  4^-)  +  1?^-  + 


3x3  v  dXi7  3Xi 


3x3  3x 


3x3'  3x5 


3x3  3xi  3Xj 


+  X,  +  A  X3  -  p  =  0  (21) 


Pi  +ipt  =o,,|i!-+p,s|Ml  +  pn  (22) 

Ps  +  4  Pi  -  on  |^  +  0,3  |Sl+  ?„  (23) 

Now  multiply  (20)  -  (23)  by  dx3  and  integrate  thru  the  initial  thickness1*; 
then  multiply  (20)  and  (22)  by  x3dx3  and  again  integrate  thru  the  thick¬ 
ness.  Into  these  results  insert  the  following  displacement  and  stress 
assumpti ons, 

Ui  (xi,  x3,  t)  =  u  (Xx,  t)  +  x3  (Xx,  t)  (24) 

u3  (xi,  x3,  t)  =  w  (xi,  t)  (25) 

an  =  E  |^-=  E  (u1  +  x 3  r)  (26) 

013  =  K2  G  (|Hl  +  |Ml)  =  K2  G  (K.  +  w')  (27) 

to  obtain  the  following  beam  equations  of  motion  and  associated  boundary 
conditions. 


t  Note  that  boundary  conditions  are  desired  only  on  the  Xx  faces,  there¬ 
fore  nx  =  1,  n2  =  n3  =  0. 

ft  that  is,  the  limits  of  all  integrals  are  from  -  h/2  to  +  h/2. 


=  p  h  u 


(28) 


(u>  Nx)#  +  (#■  Mx)'  +  (*  Qx)‘  +  E  h  u"  +  fx 
(w*  N  ) 1  +  tc2  6  h  (*'  +  w ")  +  q  =  p  h  w  (29) 

A 

CU  3 

(u'  Mx  +  i|)'  M*)'  +  (4»  Q*)'  +  "72  '  u‘  Qx  “  Q*  "^Nz  -  k2  6  h  (ip  +  W)  +  mx 

=  P  T5"  ^  (30 ) 

Fi  +  A  Fi  =  u'  N  +  \p‘  M  +  i|>  Q  +  E  h  u'  or  u  =  TT0+  (31) 

A  A  A 

F3  +  A  F3  =  w'  Nx  +  k2  G  h  (iji  +  w')  or  w  =  w0  (32) 

—  U3 

Mi  +  A  Mi  =  u'  Mx  +  3j/'  M*  +  y  Q*  +  E  V»'  or  ty  =  \p0  (33) 


where  u0,  w0,  and  tp0  have  been  used  to  denote  prescribed  displacements, 
primes  denote  differentiation  with  re/pect  to  Xi,  dots  denote  differentia¬ 
tion  with  respect  to  time,  and  where  the  following  definitions  have  been 
used. 

fx  =  /  (Xi  +  A  Xj)dx3  +  a3i^+^  (u'  +  ^-  <|>')  -  a3i^  (u '  -  |  ij>') 

+  [  a3s^  -  a33^  ]  ^  +  a3i^+^  -  o3J~^ 
q  =  /  (X3  +  A  X3) ax3  +  [  a3i^  ^  -  031^  ^  ]  w1  +  033^  ^  -  033^  ^ 
mx  =  /  (Xi  +  A  Xi)x3  dx3  +  |  [  a3i+  (u*  +  j  ij>')  +  031^  (u1  -  j  i^') 

+  (033^  +  033^)  +  a3l+  +  031^  ] 

F^  =  /  p.  dx3  (i  =1,3)  Nx  *  /  on  dx3 

A  F.  =  /  A  p.  dx3  (i  =1,3)  Nz  =  /  033 dx3 

t  Note  that  there  is  a  choice  of  stress  or  displacement  boundary  conditions. 


8 


Hi  =  /  pi  x3  dx3 
A  Mi  =  /  A  px  x3  dx3 


Qx  =  /  ffi3  dx3 


Q*  =  /  0,1  x3  dx3 


M  =  /  on  x3  dx3 


M*  =  f  o u  x3  dx 3 


Note  that  the  notation  a33^+  =  a33  (h/2)  and  o33^  ^  o33  (-  h/2)  etc.  In 


the  following  sections  the  above  results  are  used  to  analyze  some  stability 
and  vibration  problems. 


STABILITY  PROBLEMS 


This  section  employs  static  methods  (i.e.  the  "Euler  Method")  to 
determine  the  stability  of  the  various  beam  configurations  considered.  This 
is  an  adequate  approach  if  the  loadings  are  conservative.  Upon  careful 
examination  it  can  be  shown  that  the  loading  in  Case  III  is  non-conserva¬ 
tive,  hence  the  confirmation  of  that  result  must  be  postponed  until  the 
corresponding  dynamic  formulation  of  that  problem  is  carried  out  in  the 
next  section  entitled  "VIBRATION  PROBLEMS". 


Case  I,  Just  On  =  Acting 


ft 


Equations  (28)  thru  (33)  reduce  to 


(aN  +  E)  u"  =  0 

(oN  +  k2  G)  w"  +  k2  G  =  0 

^aN  +  -  k2  G  (ip  +  w' )  =  0 

u  (y)  =  ty'  (y)  =  w  (y)  =  0  ;  y  -  0,2. 


(34) 

(35) 

(36) 

(37) 


t  is  a  constant 


ft  In  all  following  examples,  the  boundary  conditions  are  chosen  to  yield 
the  simplest  solutions. 


1 

1 

§ 


I 

Js 


>3 


$ 

i 


% 


I 


pa 


I 


jV^T^V  "^s 


Neglecting  (34),  which  leads  to  trivial  results,  it  is  assumed  that 


m 


w  =  W  sin  it  x/  £ 


ip  -  ¥  cos  n  x/  £ 


whence  a  non-trivial  solution  for  W  and  ¥  demands  that  dec  | a  -  - 1  =  0 

*  J 

where 


is  I 


an  =  (S  a  -  1)  u/£ 


?i  2  =  -  1 


a2i  =  -  ir/£ 


a22  =  (o  -  K)  (S/K)  -  1 


S  =  (E/G)  (h/£)2  K  =  (£/<h)2  k2  =  tt2/12  a  =  -  aN  K/E  =  -7^-4- 

"  {  VEuler 

Hence  the  following  equation  for  a  is  obtained  where  the  lowest  root  is 


desired 


So2  -  (SK+K-f-l)a  +  K  =  0 


When  S=0  (i.e.  G  -*■  ®)  the  result  reduces  to  o0  =  K/(K  +  1).  Since  usually 


K  >  100  this  result  is  seen  to  be  close  to  the  result  obtained  from  a 
Bemoulli-Euler  beam  analysis  (o0  =  1).  For  S  ^0,  the  non-dimensional 
buckling  stress  a  is  shown  in  Figure  3  where  it  is  compared  with  the 
classical  result  *  obtained  in  Reference  5.  The  close  agreement  justifies 
the  use  of  the  simpler  analysis  developed  in  Reference  5. 


Case  II,  dust  an  =  2  x3  aM/h 


Equations  (28)  thru  (33)  reduce  to 


E  h  u"  +  M  V  =  0 


k2  G  h  (r  +  wn)  =  0 


t  The  classical  result  is  o  =  (1  +  S)" 
ft  a^  is  a  constant 


5SS Stea-feasfc* 


g 


s 


•■V  '  •~$-;;',«V^ 


<i>"  +  M  u"  -  K2  G  h  (ij»  +  w*)  =  0 


(41) 


My  *J>‘  (y)  +  E  h  u‘  (y)  =  M  u‘  (y)  +  ^  i|>'  (y)  =  w  (y)  =  0  ;  y  =  0,  l  (42) 


12 


where  Mx  =  2  h2  a^/12.  Assuming  that. 


w  =  W  sin  it  x/  £ 

'4>  =  ¥  COS  ir  x/  2. 

u  =  U  cos  it  x/  £ 


a  non-trivial  solution  for  W,  V,  and  U  demands  that  det  ja..|  =  0  where 

*  J 


an  -Eh 


a2i  -  0 


3:2  -  0 

a22  =  k2  6  h  (it/£) 


a  i  3  =  M„ 


a3i  =  Hx  (ir/£)2  a32  =  tc2  G  h  (tt/£.) 


323  -  k2  G  h 

a33  =  ic2Gh  +  Eh3  (tt/£)2/12 


Expanding  det  |a--|  =  0  yields  the  following  result, 

*  J 


°M  =  ± 


^  E 


(43) 


Thus  a  buckling  stress  is  predicted,  even  though  it  occurs  at  a  stress 
which  is  larger  than  that  permissible  for  linear  elasticity.  However,  a 
simple  plastic  yield  stress  may  be  obtained  by  substituting  the  tangent 


modulus  Et  for  E  in  (43) 


Case  III,  Just  ai3  =  Acting 


Equations  (28)  thru  (33)  reduce  to 


2  Oq  ip1  +  E  u"  =  0 


k2  G  h  +  w")  =  0 


(44) 

(45) 


t  Oq  is  a  constant 
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(46) 


££  ip"  -  k2  G  (ip  +  w')  =  0 

u  (y)  =  <P'  (y)  =  w  (y)  =  0  ;  y  =  CM  (47) 

Assuming  that, 

u  =  U  sin  it  x/  £ 

Ip  =  ¥  COS  TT  x/  £ 
w  =  Wsin-jrx/£ 

a  non-trivial  solution  for  W,  V,  and  U  demands  that  det  | a  -  - 1  =  0  where 

*  0 

3 1 1  =  E(it/£)  ai2  =  0  ai 3  =  2  Oq 

a  2 1  =  0  322  =  ^/£  323  =  1 

a3i  =0  a32  =  k2  G  (t r/£)  a33  =  <2  G  +  E  h2  (tt/£)2/12 

Expanding  the  det  |a..|  yields 

*  J 

det  la- -|  =  E2  h2  x2  (tr/£)6  G/12  j*  0  (48) 

*  J 

hence  (48)  demonstrates  that  a  buckling  stress  does  not  exist. ^  The  next 
case  demonstrates,  however,  that  if  Oq  is  applied  conservatively,  a  buck¬ 
ling  stress  does  exist. 

Case  IV,  Just  aJ3  =  aQ  Acting  Conservatively 

In  this  case  f  -  m  =  0  so  that  (28)  thru  (33)  reduce  to 

X  X 

aQ  ip'  +  E  u"  =  0  -  (49) 


t  The  confirmation  of  this  result  is  found  in  Case  IV  of  the  next  section 
entitled  "VIBRATION  PROBLEMS". 
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ic2  G  h  (tfi‘  *  w")  =  0 


E  V;"  -  0^  u‘  -  tc2  G  (^  +  w‘)  =  0 


u  (y)  =  <J>‘  (y)  =  w  (y)  -  0  ;  y  -  0,5. 


Assuming  that 


u  =  U  sin  tt  x/  5 


=  ¥  cos  it  x/  5 


w  =  W  sin  it  x/  £ 


a  noti-trivial  solution  for  W,  'F,  and  U  demands  that  det  |a..|  =  0  where 

J 


an  =  K  (tt/£) 


a  1 2  =  0 


ai 3  =  K  Oq/E 


a2i  ~  0 


a2-i  =  (it/£)/S  a23  =  1/S 


a3i  =  K  Oq  (tt/£)/E  a32  =  (it/£)/S  a33  =  1  +  1/S 


Expanding  det  j  a  -  - 1  =  0  yields  the  following  result, 

•  J 


Oq  =  ±  E  K  (h/£) 


This  result  has  been  determined  previously  by  Herrmann  and  Armenakas  [2]. 


However,  as  demonstrated  by  Case  III,  their  statement  "...  we  could  assume 


that  the  resulting  shear  forces  either  rotate  with  the  element  or  do  not 


change  direction  after  deformation."  can  not  be  correct. 


VIBRATION  PROBLEMS 


Case  I,  Longitudinal  Vibration  with  Just  an  =  aM  Acting 


In  this  case  (28)  and  (31)  are  uncoupled  and  become 


(oN  +  E)  u"  =  p  il 

u  (y,t)  =0  ;  y  =  0,£ 

Assuming  u  =  Un  (sin  mr  x/  i)  elwn*  and  letting 
is  given  by 


(54) 

(55) 

±  a  E  the  solution 


wn 

Ti-7 


Cl 


(1  ± 


(55) 


where 

(to  )  =  (E/p)1^2  n  7r/  & 

n  cl 

and  1 

n  =  1 ,  2,  3,  ... 

For  some  materials  (glass  fibers,  for  example)  a  may  be  as  larjje  as  .1, 
hence  for  this  value  of  a  it  is  seen  that  (56)  may  be  in  the  range  from 
.95  to  1.05  for  all  values  of  n. 


Cas~  II,  Transverse  Vibration  with  Just  an  =  a^  Acting 


In  this  case  just  (29),  (30)  and  (32),  (33)  are  coupled  and  become 


(a^  +  k2  G)  w"  +  k2  G  iff'  =  p  w 


(aN  +  E)  ip"  -  k2  G  (ij>  +  w')  =  p  Yp  \p 


2  - 
12 


ij>'  (y,t)  =  w  (y,t)  =0  ;  y  =  0,£ 


Assuming  that. 


w  =  Wn  (sin  n  x  x/  £)  ela3nt 


=  >},n  (cos  n  i  x/  £)  e1wn* 
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(57) 


(58) 

(59) 


a  non-trivial  solution  for  Wp  and  ^  demands  that  det  [a^l  =0  where 

an  =  n2  (S  o  -  1)  +  fl2n  S  ai2  =  -  n  JI/tt 

a2i  =  -  n  K  r/£  a22  =  n2  S  (a  -  K)  -  K  +  S  ft2n 

and 

fl2n  =  p  K  (A/it)2  a)2n/E 
The  frequency  spectrum  thus  is  given  by 

fln  =  [  (-  B  ±  [B2  -  4  A  C]1/2  )/2  A  ]1/2  (60) 

where 

A  =  S 

B  -•  n2  o  S  +  n2  S  (a  -  K)  -  K  -  n2 
C  =  n2  a  [  n2  S  (a  -  K)  -  K]  -  n4  (a  -  K) 

When  S  -*•  0  (zero  transverse  shear)  (60)  specializes  to 

=  f  (-  n2  a  K  +  n4  [K  -  a]  )/[K  +  n2]  ]1/2  (61) 

Note  that  when  a  =  S  =  0  and  K  is  large,  fin  -  n2  which  is  the  classical 
result.  Figure  4  plots  the  first  two  non-dimensional  frequencies  versus 
a  for  three  values  of  E/G^  with  £/h  =  10.  The  classical  values  [5]  are 
shown  by  circles  and  are  in  very  good  agreement  with  the  present  results, 
thus  completely  confirming  the  results  obtained  in  Reference  5. 

Case  III,  Vibration  with  Just  an  =  2  x3  a^/h  Acting 
Equations  (28)  thru  (33)  reduce  to 


t  Note  that  E/G  =  0  yields  the  classical  values  and  that  E/G  =2.6  yields 
the  isotropic  values. 


IS 


M  y  +  E  h  u"  =  p  h  u 

A 


k2  G  h  (ty1  +  w")  =  p  h  w 


Mx  u“  +  ^*u  -  k2  G  h  (if.  +  w1)  =  p  !£* 


(62) 

(63) 

(64) 


Mx  (y,t)  +  E  h  u'  (y,t)  =  My  u'  (y,t)  +  ^  V  (y,t) 


=  w  (y,t)  =0  ;  y  =  0,£  (65) 


Assuming  that 


w  =  Wn  (sin  n  it  x/  i)  ela)n^ 

ip  =  ^  (cos  n  tr  x/  £)  e1tl3n^ 

u  =  Un  (cos  n  it  x/  £)  e1a)n* 


a  non- trivial  solution  for  Wn>  Yn,  and  Un  demands  that  det  la^j  =0  where 


a  =  K  Oft/ E  and 


an  =  K  n2  -  fi2n  an  -  0  a13  =  2  n2  £2  d/K  h  tt2 


a2i  -  0 


a22  =  nz  -  ft2  S  a23  =  n  £/tt 


a3i  =  2  n2  S  o/h  a32  =  n  it  K/£  a33  =  K  (1  +  n2  S)  -  S  ft2 


Expansion  of  the  determinant  yields  the  following  expression  for  ft 


(K  n2  -  ft2n)  [  n1*  K  -  ft2n  (n2  +  K  +  n2  K  S)  -  S  ft^  ]  +  (ft2n  S  -  n2)  (n2  o)2/3  =  0  (66) 

|  For  zero  transverse  shear  S  -*■  0  and  (66)  reduces  to 

! 

(n2  +  K)  ft"n  -  K  n2  (2  n2  +  K)  ft2p  +  n6  (K2  -  a2/3)  =  0  (67) 

Noting  that  a  appears  only  in  the  zero^  power  term  in  (66)  and  (67)  and 
noting  that  <  .1  E  (hence  a  <  .1  K)  it  is  seen  that  a  has  a  negligible 
effect  on  the  frequencies. 
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Case  IV,  Vibration  with  Just  g13  =  Acting 


Equations  (28)  thru  (33)  reduce  to 


2  Oq  4>‘  +  E  u‘ 


■"  -  p  u 
K2  6  (l|>‘  +  w")  =  P  V/ 


E^r  -  .c2  6  (*  +  w*)  =  p^* 


U  (y,t)  =  If/’  (y,t)  =  w  (y,t)  =  0  ;  y  =  0,2, 


(68) 

(69) 

(70) 

(71) 


Assuming  that 


u  =  Un  (sin  n  tt  x/  £)  elton^ 


iJj  =  ¥  (cos  n  ir  x/  £)  ela3n^ 


w  =  Wn  (sin  n  n  x/  £)  e1U)n^ 


a  non-trivial  solution  for  U  ,  ¥  ,  and  W  demands  that  det  ja..|  =  0  where 

n  n  n  1  lj1 

a*  =  -  K  an/E  and 


Q 


au  =  fl2n  -  K  n2 


di2  -  0 


au  =  -  2  o*  n  £/tt 


a2i  -  0 
a3i  ~  0 


a22  =  S  fi2n  -  n2  a23  =  -  n  l/tr 


a32  =  -  n  it  K/Z  a33  =  S  ft2n  -  K  (1  +  n2  S) 


Expansion  of  the  determinant,  given  below,  demonstrates  that  a*  does  not 


affect  any  of  the  frequencies  and  in  fact  the  longitudinal  frequency  is 
uncoupled  from  the  other  two  coupled  frequencies. 


(«2n  -  K  n2)  [  (S  fi2n  -  n2)  (S  n2n  -  K  [1  +  n2  S])  -  n2  K  ]  =  0  (72) 


For  zero  transverse  shear  S  -*■  0  and  (72)  reduces  to 


t  In  particular  no  value  of  a*  can  make  the  frequencies  zero  or  complex, 
hence  instability  cannot  occur  as  was  predicted  in  Case  III  of  the  pre¬ 
vious  section  entitled  "STABILITY  PROBLEMS". 


17 


js 

% 


■a 


% 


45 


i 


3 

1 


I 


I 


% 

I 


J 


(79) 


%=  [  [-  B  ±  <B2  -  4  A  C)1/2]/2  A  ]1/2 

where 

A  =  K  +  n2 

B  =  -  [  n4  K  +  n2  K  (K  +  n2)  ] 

C  =  n4  K  [  n2  K  -  a2  ] 

Figure  5  plots  the  first  two  non-dimensional  frequencies  versus  a*  for 
three  values  of  E/6  with  £/h  -  10. 

COMPARISON  WITH  A  PREVIOUS  THEORY 

f 

Referring  to  Reference  2  and  one-dimensionalizing  the  results  yields 
the  followir.g  equations  to  be  compared  with  (28),  (29),  and  (30). 

(¥  Qx)'  +  E  h  u"  +  fx  =  p  h  u  (80) 

(w1  Nx)'  +  k2  G  h  (f  +  w")  +  q"  =  p  h  w  (81) 

■^2  f  0x-K2Gh(t  +  w')+inx  =  p^  (82) 

It  is  noticed  that  the  fonn  of  (81)  and  (29)  agree,  but  that  the  forms  of 

(80)  and  (28)  and  of  (82)  and  (30)  have  marked  differences.  Additionally, 
Reference  2  leaves  the  terms  7  ,  q”,  and  mx  to  be  arbitrary  whereas  the 
present  paper  demands  specific  forms  for  f  ,  q,  and  m  as  seen  by  the  first 

three  results  following  (33).  Also  in  a  previous  criticism  Masur  [6]  has 

pointed  out  that  a  term  ( u '  Nx)'  should  appear  in  (80)  which  is  accounted 
for  in  the  corresponding  equation  (28)  of  the  present  work.  Finally,  since 


t  Notational  changes  have  been  made  so  that  the  comparison  of  equations  (36a), 
(36c),  and  (36d)  of  Reference  2  with  (28),  (29),  and  (30)  of  the  present 
work  may  be  facilitated. 


I 


m 


the  present  work  is  based  on  simple  straight-forward  origins  and  is  so 
easy  to  check  it  is  difficult  to  imagine  that  it  can  be  in  error;  hence 
one  is  forced  to  the  conclusion  that  the  one-^imensionalized  results  of 
Reference  2  are  not  correct. 


CONCLUSIONS 

The  several  conclusions  that  may  be  inferred  from  this  investigation 
are  that  (a)  with  just  acting,  the  simpler  work  of  Reference  5  is  ade¬ 
quate  to  describe  the  problem,  (b)  with  just  <Jq  acting,  very  small  changes 
in  the  direction  of  loading  exert  a  large  influence  on  determining  whether 
or  not  the  system  can  be  unstable,  (c)  with  just  aM  acting  negligible  ef¬ 
fects  are  found,  (d)  transverse  isotropy  still  exerts  a  significant  influ¬ 
ence  on  stability  and  vibration  regardless  of  the  type  of  initial  stresses 
that  are  acting,  and  finally  (e)  the  one-dimensionalized  version  of  a  pre¬ 
vious  work  [2]  by  Herrmann  and  Armenakas  does  not  accurately  describe  the 
present  problem,  hence  doubt  is  cast  as  to  the  validity  of  their  paper  in 
general . 
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NOMENCLATURE 


&£ 

R,* 

Sr- 


1 

I 

1 


& 


I 


E 

^  +  A  Fi  i  i  »  1,3 

V  <*»  mx 
G 

h 

K  =  (  JL/kH  ) 2 
l 

Mi  +  A  Mi 

Nx,  Nz,  Qx.  Q*.  Mx,  M* 
P* 

PS 

P-j  +Api 

S  -  (E/G)  (h/£)2 
u,  w 
TTo,  w0 


longitudinal  Young's  modulus 
perturbation  boundary  forces 
perturbation  loadings 
transverse  shear  modulus 
beam  thickness 
geometry  parameter 
beam  length 

perturbation  boundary  moment 
initial  generalized  stress  resultants 
generalized  traction  vector 
actual  traction  components 
perturbation  boundary  stresses 
transverse  isotropy  parameter 
perturbation  displacements 
perturbation  displacements  at  the  boundary 


generalized  body  force  vector 
actual  body  force  components 
final  quantities 
perturbation  quantities 
Mindlin's  shear  correction  factor 
perturbation  rotation 
perturbation  rotation  at  the  boundary 
generalized  stress  vector 
perturbation  stresses 
non-dimensional ized  shear  stress 
non-dimensional ized  outer  fiber  stress 


non-dimensionalized  normal  stress 
constant  initial  stresses 
Trefftz  components  of  stress 
actual  stresses 
non-dimensional  frequency 
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